Single hole doped strongly correlated ladder with a static 

impurity 
S. Gayen 

Department of Physics, Drexel University, Philadelphia, PA 19104, USA 



Abstract 



We consider a strongly correlated ladder with diagonal hopping and exchange 
interactions described by t — J type hamiltonian. We study the dynamics 
of a single hole in this model in the presence of a static non-magnetic (or 
magnetic) impurity. In the case of a non-magnetic (NM) impurity we solve 
the problem analytically both in the triplet (S=l) and singlet (S=0) sectors. 
In the triplet sector the hole doesn't form any bound state with the impu- 
rity. However, in the singlet sector the hole forms bound states of different 
symmetries with increasing J/t values. Binding energies of those impurity- 
hole bound states are compared with the binding energy of a pair of holes 
in absence of any impurity. In the case of magnetic impurity the analytical 
eigenvalue equations are solved for a large (50 X 2) lattice. In this case also, 
with increasing J/t values, impurity- hole bound states of different symme- 
tries are obtained. Binding of the hole with the impurity is favoured for 
the case of a ferromagnetic (FM) impurity than in the case of antiferromag- 
netic (AFM) impurity. However binding energy is found to be maximum for 
the NM impurity. Comparison of binding energies and various impurity-hole 
correlation functions indicates a pair breaking mechanism by NM impurity. 



1 



1. Introduction 



Since the discovery of the high-T c superconductivity [T] the nature of the 
order parameter, the symmetry of the bound pair of holes has been a hotly 
debated issue [2-10]. It is well known that in the case of conventional su- 
perconductors the introduction of magnetic impurities drastically suppresses 
the transition temperature T c . Also theoretically it was predicted that mag- 
netic impurities are strong pair breakers for conventional superconductors 
while NM impurities have a pair breaking effect only for higher orbital 
momentum states [T2*l IT3] such as d-wave pairing state. There were exten- 
sive studies on the effect of doping the high T c cuprate systems with static 
magnetic and NM impurities in the last decade [14-17]. Experimental stud- 
ies have shown that divalent Zn and Ni ions replace the Cu ++ ions in the 
CuC>2 planes and ZN (the NM impurity) is more effective than Ni (the mag- 
netic impurity) in destroying superconductivity. It was suggested by some 
groups E3 that NM impurities do not simply act as vacancies but in- 

duce local magnetic moments in the CuC>2 planes and the mechanism is still a 
magnetic pair breaking mechanism. However, other experimental studies [2U] 
have indicated that the estimated impurity moment-carrier exchange in Zn 
doped YBCO is too small to account for the suppression of T c by a magnetic 
pair breaking mechanism satisfactorily. Moreover, the observed difference in 
behaviour of Zn and Ni doped systems could not be properly explained. 

In order to explain all these features qualitatively Poilblanc etal. initiated 
numerical studies on a microscopic model, the t — J model. In a series of 
papers [2U 122 E3] they studied the effect of magnetic and NM impurities 
in very small clusters in 2D. Based upon their calculation of local density 
of states in presence of an impurity (NM or magnetic), they found out the 
existence of successive impurity-hole bound states of d,s and p-wave symme- 
tries with increasing J/t values [2BI22]- From the analysis of impurity-hole 
binding energies and quasiparticle weight of a pair of holes [2B] they con- 
cluded that NM impurity has a stronger pair breaking effect than a magnetic 
impurity. In this respect we also looked at the effect of magnetic and NM 
impurities on single hole dynamics in a t — J ladder model introduced by 
Bose and Gayen [21]. Existence of bound hole pair and superconductiong 
correlations in a t — J ladder model was first numericallly shown by Dagotto 
etal. in 1992 [23]. Later Gayen and Bose j2Hl I2Z| proved analytically that 
the two holes form a bound propagating object in the ground state with a 
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modified d-wave symmetry and the superconducting correlations exist in a 
t — J ladder model with diagonal interactions. 

We divide the paper in the following sections. In section 2 we briefly 
introduce the model and recapitulate the exact results in the undoped limit 
as well as for single and two holes doped cases. In section 3 we derive 
analytical results for single hole dynamics in presence of a static NM impurity. 
In section 4 we consider the hole in presence of a static magnetic impurity. 
In the last section (section 5) we conclude by speculating the implication of 
a pair-breaking mechanism by NM impurities in our model. 

2. The Model 

The ladder model shown in figure 1 consists of two chains joined together by 
rungs and diagonals. The model is described by the following hamiltonian 

h = - E '■M'.J'j^ ■ "■('■) ■ E JiA-Si (i) 

<i,j>a <i,j> 

The hopping integral along the solid (dashed) line is t (£') . The corresponding 
AFM exchange interaction is J (J'). The spins are of magnitude | and 
single occupancy constraint is strictly enforced in our model. We use periodic 
boundary conditions (PBC). In the undoped limit the ground state consists of 
singlets along the vertical rungs when 2 J = J' [21] with ground state energy 
E G = -\J'N. Actually this is the ground state for J' > 1.40148403897 J 
[28J. In this paper we shall consider J' = 2 J for simplicity. However, the 
results are exact for all values of J' > 1.40148403897J. When the ladder is 
doped by a single hole, depending upon the value of J/t the ground state 
may consist of either the hole bound to an adjacent triplet excitation or the 
hole can coherently propagate through a sea of singlets [2H]. When the ladder 
is doped with two holes the ground state consists of coherently propagating 
bound state of two holes of modified d-wave symmetry [2"Ej . 

In this context it is important to note that ladder models are not only of 
theoretical interest but also of experimental relevance [3U]. Hiroi etal. [3*T] 
synthesized the family of compounds 5'r n _ 1 Cu n+1 2ri that can be described 
well by ladder geometry [321 EH] • Midler and Mikeska [33j discussed about 
the elementary excitations in the compound KCuCl% which is described by 
an AFM ladder model with additional diagonal couplings. The compound 
LaCu0 2 .5, formed by an array of weakly interacting ladders [US], could be 
doped with holes [HH] on replacing La by Sr. Furthermore, the discovery of 
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superconductivity in doped ladder compound Sr\^ x C a x Cu2AO ^\ by Uehara 
etal. [H7| gave a tremendous impetus to theoretical and experimental research 
works on ladders. 



Figure 1: The ladder model with diagonal coupling. 

3. Non-magnetic Impurity 

The NM impurity is modelled as an inert site located on the first rung of 
the ladder. PBC is implied. The ground state consists of singlets along all 
the rungs except the first rung where there is a free spin-| along with the 
impurity. For N rungs the ground state energy will be E g = —^j-(N — 1). 
A single hole can be added in the singlet state by removing the free spin. It 
can be added to another rung thereby creating another free spin also. The 
two free spins in the first rung and in the other rung, where a hole had been 
added, can combine to form either a singlet (S=0) or triplet (S=l) state. 

3.1 triplet (S=l) sector 

In the triplet sector we can generate a closed subspace of basis states 
0(l),0(2),---,0(iV-l), 

0(1) = t-4(J)+?)ll-*' 



0(2) = f|-^(2>+?) 



V2' 
1 

1 

3 V 



,ir-s(J>+?)ll- (2) 



where the argument r implies that the hole is on the (r + l)-th rung. The ver- 
tical solid line along i-th rung stands for a singlet (^(c 2 j_ijc 2 ^ — c 2 j_i^c 2i f )). 
A general eigenfunction \1> = Y.fji 1 a i<K0 with energy eigenvalue E satisfes 
the following eigenvalue equations: 



—)a x = ta 2 
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ea 2 = t(ai + a 3 ) 

eajv_ 2 = t(ciN-i + (Ins) 
(e - = ta7v_ 2 (3) 



where e = E — ^j- + t' . The energy eigenvalue i£ is measured with re- 
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spect to the ground state energy The above eigenvalue equations can be 
solved analytically with symmetric (a n = cos(k(N/2 —n))) or antisymmetric 
(a n = sin(k(N/2 — n))) eigenfunctions. In the symmetric case the energy 
eigenvalues are obtained from simultaneous solution of the equations 

e = 2tcos(k) 
_ J cos(k(N/2-2)) 
C 4 cos(k(N/2 - 1)) { ) 

In the antisymmetric case the corresponding equations are 

e = 2tcos(k) 
J sin(k(N/2 - 2)) 

C ~4 ~ sin(k(N/2 - 1)) W 

In the limit — > oo we get a continuum of extended states for real values of 
k. When J > At we can derive antibound state solutions in both the cases for 
imaginary values of k (k replaced by ik). The antibound states are localized 
states lying above the continuum of extended states with energy eigenvalue 
E a b = — ^ / + i + ~7~- However, for any positive value of J the hole does 
not form a bound state with the impurity in the triplet sector. 

3.2 singlet (S=0) sector 

The motion of hole in presence of NM impurity were reported in brief by 
Bose and Gayen [[JH] earlier. For the sake of convenience we briefly review 
the salient features here. A basis state fa will have the hole located in the 
i-th rung. There will be singlets residing on all — 2 rungs except the 
first and i-th rungs. The electronic spins on 1st and z-th rungs will form 
a singlet. In the basis state fa both the impurity and the hole are on the 
first rung. A general eigenfunction will have the form \I/ = Ylj^i^'Pi- The 
exact amplitude equations derived from the time independent Schrodinger 
equation with energy eigenvalue E (measured with respect to the ground 
state energy E g ) are 
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3 J' 

(e + — t')ai = tV2(a 2 + a N ) 



4 

e+—)a 2 = -rv2ai + ra 3 



ea 3 = t(a 2 + a 4 ) 



eaAr_i — i(a/v + ojv-2) 
3 J 

(e + — )oat = -tv^ai + taAr_i (6) 

where e = E — H- + t'. In this problem also the eigenfunctions are symmetric 
or antisymmetric with respect to reflection about the first rung containing 
the impurity. In both cases for real values of k one gets extended states. 
With appropriate choice of complex values of k one can find out bound 
(k — > 7T + ik) or antibound (k —>■ ik) states in the symmetric sector. In the 
case of antisymmetric wave function one only finds out bound state solutions. 
In the N — > 00 limit, we derive the following results: (a) For the symmetric 



eigenfunctions, (i) for J < |[|- + 2t — y8t 2 + (f-) 2 ], there can be no bound 
state solution; (ii) for J > |[~ — 2t + J8t 2 + (f-) 2 ], we do not find any 



antibound state; (iii) for J > + 2t + J8t 2 + (f-) 2 ], we always find two 
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bound states. 

(b) In the case of antisymmetric eigenfunctions, we always get one and only 
one bound state provided J > 3 ff^ 2 ) • ^ n ^ e ^ ~^ 00 nm ^> ^ ne bound state 
in the antisymmetric sector has energy E = — t' — ^- — ^j- 

We use the following quantitative definition of the impurity-hole binding 
energy introduced by Poilblanc etal. [2U 123 

A B (l/i, li) = (E(lh, li) - E(0h, li)) - Oi) - E(0h, Oi)) (7) 

Figure 2 shows the plot of binding energy (A^/t) as a function of exchange 
interaction (J/t) in different symmetry channels. We also compare the bind- 
ing energy of a pair of holes computed from the results of Gayen and Bose 
in the same figure. The binding energy of two holes in absence of impurity 
is defined as 

A B {2h, Oi) = (E(2h, Oi) - E{0h, Oi)) - 2{E{lh, Oi) - E(0h, Oi)) (8) 

where E(nh,0i) is the ground state energy in presence of n holes and 
impurity. E(lh,li) is the lowest energy in presence of one hole and one 
impurity in the particular symmetry channel. 
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O 1 2 3 4 5 6 7 8 910 



J/t 

Figure 2: Binding energy As/t vs. J/t for bound states of different symme- 
tries in the case of NM impurity. The binding energy of a pair of holes is 
also shown in the figure (t'/t = 1.0). 

From figure 2 we find that for any positive value of J / 1 there is always one 
bound state in the symmetric (symmetric 1) channel. For J/t > ~ one can 
find another bound state in the antisymmetric channel. Another bound state 
is obtained in the symmetric (symmetric 2) channel for all values of J/t > 
3.58152. The co-existence of the bound states in different symmetry channels 
are analogous to the appearance of d-, s-, p- wave bound states in 2D clusters 
studied by Poilblanc etal. [21 J . However, there is one significant difference 
in the ladder. The sequence of the bound states appearing in symmetry 
channels of different parity (even-odd-even) is different from that (even-even- 
odd) obtained by Poilblanc etal. [21]. The most significant observation from 
the graph is that the impurity-hole binding energy in the lowest symmetric 
sector has magnitude greater than the binding energy of a pair of holes for 
J/t > 0.74. For 2D finite clusters Poilblanc etal. pT[| found the binding 
energy curve for the pair of holes to be lowest. 
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Figure 3: Impurity- hole correlation Cih in the ground state as a function of 
the distance dm of the hole from the NM impurity (t = t' = 1.0). 

In order to have more understanding about the nature of impurity-hole 
binding we plot (figure 3) the impurity- hole correlation function Cih{j — 1) 
(=< rii(l)nh(j) >) as a function of the distance d^ between the impurity 
and the hole in the ground state. When the hole shares the same rung with 
the impurity, there is no extra cost of energy due to broken singlets along 
the rungs. If the hole sits one rung away from the impurity, one additional 
singlet along the rung is broken. But this loss of energy can be compensated 
by the creation of an extra singlet between the two free spins sitting on the 
two rungs which are occupied by the impurity and hole also. Moreover, the 
hole can gain kinetic energy of the amount -t' when it stays one rung away 
from the impurity. As a result the hole likes to sit one rung away from 
the impurity in the ground state for smaller values of J (with t — t' — 1). 
However, with increasing J values the probability of the hole sharing the 
same rung with the impurity increases rapidly. 

4. Magnetic Impurity 

We assume the static impurity to be a spin-| object that interacts with neigh- 
bouring spins via exchange interaction J' m (along the rung) and J m (along 
horizontal and diagonal bonds). The dynamics of the system is governed by 
Himp which is essentially the hamiltonian introduced in eqn.(l) with the fact 
that the electrons can not hop on to the impurity site and the spin at the im- 
purity site interacts with the neighbouring spins with different strengths. We 
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assume that the static impurity sits on the first rung. We use PBC. In the ab- 
sence of any hole the ground state consists of singlets along all the rungs with 
energy eigenvalue —(^j-(N — 1) + ^f^) if the impurity exchange interaction 
is AFM. For FM exchange interaction, all the rungs, except the first one, are 
singlets. The first rung, where the magnetic impurity is located, has 5 = 1 
configuration. The corresponding energy eigenvalue is —(^j-(N — 1) + -^p). 
A single hole can be doped in the first rung to create the state 

X=SH---- 0) 

We can diagonalize the hamiltonian within a closed subspace of basis states 
denoted by, X , 0(1), 0(2), • • •, </>(N-l), 9(1), 9(2), ■ ■ ■, 9(N - 1), ^(1), ^(2), 
• • •, ip(N — 1). The states are pictorially represented in the following way 



0(r) = 




■ ■ (?+ l o) 


9(r) = 


(^)rn- 


■ ■ (?+'o) 


ip(r) = 


(>n- 


■ ■ (?+h) 



where the argument r implies that the hole is located 
A general eigenfunction will have the following form 



(10) 

in the (r + l)-th rung. 



N-l N-l N-l 

* = a oX + E b r4(r) + £ c r 9(r) + £ d^(r) (11) 

r=l r=l r=l 

Clearly the total number of basis states will be mm = 1+3(N —1) and we have 
to diagonalize mm x mm hamiltonian to find out the eigenvalues and eigen- 
functions. In our numerical calculations we take J' m = 2J m , in accordance 
with the interactions (J' = 2 J) in between the host spins. In this case also the 
eigenfunctions are found to be symmetric (b n = biy-n, Cn = Qv-n, d n = djsr-n) 
or antisymmetric (b n = —bN^ n ,c n = —c^- n ,d n = — djv-n)- The lowest en- 
ergy eigenstate is always obtained in the symmetric sector for both FM ar 
AFM impurity. In figure 4 we compare the binding energy of bound states 
in different symmetry channels for FM, AFM and NM impurities. The re- 
sults obtained indicate that the binding of a doped hole with the impurity is 
favoured in the case of FM impurity. Binding occurs for weakly AFM impu- 
rity also. However, the binding energy in a particular symmetry channel is 
maximum if it is a NM impurity. In figure 5 we plot the binding energy in 
the ground states for FM (J m = —0.5) and AFM (J m = 0.5) impurity and 
compare them with that in the case of NM impurity. 
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symmetric 1 
antisymmetric 
symmetric 2 



Figure 4: Binding energy of bound states of different symmetries as a function 
of impurity exchange coupling J m (J = 5.0,£ = t' = 1.0, N = 50). 




j 

Figure 5: Binding anergy of the ground state as a function of J/t for FM 
(J m = -0.5), AFM (J m = 0.5) and NM (J m = 0.0) impurity (* = *' = 1.0, 
iV = 50). 

In figure 6 we plot the impurity-hole correlation function (Cih) for dif- 
ferent types of impurities. We notice that the impurity-hole correlation is 
maximum when the hole is just one rung away from the impurity. This is ex- 
pected from figure 3 for the present choice of parameters (J = t = t') . When 
the hole is one rung away from the impurity, all other parameters remaining 
the same, the magnitude of the correlation function is maximum for AFM 
impurity and minimum for FM impurity. Actually For AFM impurity all the 
three basis states (f)(1), 9(1) and ip(l) (all of them have the hole occupying 
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the first rung) contributes to a large extent to the correlation function. For 
FM impurity there is little contribution form the basis state ^(1) due to un- 
favourable configuration (two of the three bonds are frustrated). However, 
if the hole and impurity are on the same rung, the correlation function has 
largest value for NM impurity and smallest value for AFM impurity. 




d , h 



Figure 6: Impurity- hole correlation Cih in the ground state as a function of 
the distance dih of the hole from the impurity (FM, AFM and NM) (J = t = 
i! = 1.0, N = 50). 

5. Conclusions 

In this paper we have solved analytically the many-body problem of a single 
hole in presence of a NM impurity and infinite number of electrons for a 
ladder of infinite length. A continuum of scattering states and antibound 
state are obtained in the triplet sector. In the singlet sector, in addition to the 
states found in the triplet sector, bound states of the hole with the impurity 
are also found. As the ratio J/t is increased from zero, successive bound 
states in the symmetric, antisymmetric and symmetric channel appear. For 
large values of J/t all the three bound states co-exist. As pointed out by 
Poilblanc etal. 21\ this is a characteristic feature of strong correlation. Our 
analytical calculations for the ladder model can qualitatively confirm the 
numerical results for small 2D clusters by Poilblanc etal. |21[ I22j. However, 
contrary to the observation by Poilblanc etal. we find that the binding energy 
curve for the impurity-hole bound state (in the symmetric sector) is lower 
than that of a bound pair of holes in almost the whole of the parameter 
space and this definitely indicates a strong pair breaking effect [39J by NM 
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impurity in the t — J ladder model. In the case of single hole in presence 
of a magnetic impurity, the eigenvalue problem is solved numerically for a 
finite-sized ladder. Since the number of basis states increases as N we can 
solve a large size ladder with relative ease. Again, the results obtained are 
in qualitative agreement with those obtained by Poilblanc etal. [22] for 2D 
finite clusters. Impurity-hole binding is more favourable in the case of FM 
impurity than in the case of AFM inpurity. With increase in the strength of 
the AFM interaction of the impurity with surrounding electronic spins, the 
magnitude of binding energy rapidly decreases to zero. Comparison of the 
results between magnetic and NM impurites indicates that NM impurities 
should act as stronger pair-breaking agent in strongly correlated systems 
than magnetic impurities. 
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